POLYNOMIAL APPROXIMATION IN SOBOLEV SPACES ON 
THE UNIT SPHERE AND THE UNIT BALL 



FENG DAI AND YUAN XU 

Abstract. This work is a continuation of the recent study by the authors on 
approximation theory over the sphere and the ball. The main results define 
new Sobolev spaces on these domains and study polynomial approximations 
for functions in these spaces, including simultaneous approximation by poly- 
nomials and relation between best approximation to a function and to its 
derivatives. 



1. Introduction 

In a recent work [3], the authors defined new moduli of smoothness and K- 
functionals on the unit sphere — {x G R d : ||x|| = 1} and the unit ball 

M d = {x G R d : ||x|| < 1} of R d , where ||x|| denotes the usual Euclidean norm, and 
used them to characterize the best approximation by polynomials. This work is a 
continuation of [3] and studies polynomial approximation in Sobolev spaces. 

The new modulus of smoothness on the sphere is defined in terms of forward 
differences in the angle 6ij of the polar coordinates on the (x{,Xj) planes, and it 
is essentially the maximum over all possible angles of the moduli of smoothness 
of one variable in O^j. There are ( 2 ) such angles, which are clearly redundant as 
a coordinate system. Nevertheless, our new definition effectively reduces a large 
part of problems in approximation theory on S d to problems on S 1 , which allows us 
to tap into the rich resources of trigonometric approximation theory for ideas and 
tools and adopt them for problems on the sphere. These same angles also become 
indispensable for our new definition of moduli of smoothness on the unit ball M d . 
In fact, our moduli on M d are defined as the maximum of moduli of smoothness 
of one variable in these angles and of one additional term that takes care of the 
boundary behavior. We had two ways to define the additional term, the first one is 
deduced from the results on the sphere and the second one is the direct extension 
of the Ditzian-Totik modulus of smoothness on [—1,1], both of which capture the 
boundary behavior of the unit ball and permit both direct and inverse theorem for 
the best approximation. For d = 1, approximation by polynomials on B 1 = [—1, 1] 
is often deduced from approximation by trigonometric polynomials on the circle 
S 1 by projecting even functions and their approximations on S 1 to [—1,1]. This 
procedure can be adopted to higher dimension by projection functions on E> d onto 
M d , and this is how our first modulus of smoothness on M d was defined. It should 
be mentioned that our new moduli of smoothness on the sphere and on the ball are 
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computable; in fact, the computation is not much harder than what is needed for 
computing classical modulus of smoothness of one variable. A number of examples 
were given in [3J. 

In the present paper we continue the work in this direction, study best approx- 
imation of functions and their derivatives. On the sphere, our result will be given 
in terms of differential operators Dij = x^dj — Xjdi, which can be identified as par- 
tial derivatives with respect to 6i.j. We shall define Sobolev spaces and Lipschitz 
spaces in terms of D^j on these two domains and study approximation by polyno- 
mials in these spaces, including simultaneous approximation of functions and their 
derivatives. The study is motivated by a question of Kendall Atkinson (cf. [1]) in 
numerical solution of Poisson equation. 

The paper is organized as follows. The main results in j3j on the unit sphere 
will be recalled in Section 2 and the new results on the sphere will be developed in 
Section 3. The results in [3J on the unit ball will be recollected and further clarified 
in Section 4. Finally, the new results on the ball are developed in Section 5. 

Throughout this paper we denote by c, Ci,C2,... generic constants that may de- 
pend on fixed parameters and their values may vary from line to line. We write 
A < B if A < cB and A ~ B if A < B and B < A. 



2. Polynomial approximation on the sphere: Recent Progress 

In this section we recall recent progress on polynomial approximation on the 
sphere as developed in [3J. 

Let L p (S d_1 ) be the L p -space with respect to the usual Lebesgue measure da 
on with norm denoted by || ■ || p := || • ||lp(§<i-i) for 1 < p < oo, where for 

p = oo, we replace L°° by C(S d_1 ), the space of continuous functions on E> d ~ 1 with 
the uniform norm. 

2.1. Polynomial spaces and spherical harmonics. We denote by LT^ the space 
of polynomials of total degree n in d variables, and by II n (S *) := II„|gd-i the 
space of all polynomials in II„ restricted on In the following we shall write 11^ 

for nf^S^ -1 ) whenever it causes no confusing. The quantity of best approximation 
is then defined by 

(2.1) E n (f) p := inf ||/- 9 ||p, 1 < p < oo. 

Let T-V^ denote the space of spherical harmonics of degree n on S d_1 , which are 
the restrictions of homogeneous harmonic polynomials to § d_1 . Let Ao be the 
Laplacc-Beltrami operator on the sphere, defined by 



(2.2) A /(a;) 



/ 



V 

h\ 



d-l 



where A := + . . . + is the usual Laplace operator and it acts on the variables 
y. Then the spherical harmonics are the eigenfunctions of Ao, 

(2.3) A r = ~n{n + d - 2)Y, Y eUt 

The reproducing kernel of the space T-L^ in L 2 (E> d ~ 1 ) is given by the zonal harmonic 

(2.4) Z n<d (x,y) := ^±±C*((x,y)), A = 
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where (x, y) denotes the Euclidean dot product of x, y £ R d . Let be the Gegen- 
bauer polynomial with index A, normalized by C*(l) = ~ )• Let 77 be a C°°- 
function on [0, 00) with the properties that rj(x) = 1 for < x < 1 and r)(x) = 
for x > 2. We define 



(2.5) V n f(x):= f(y)K n ((x >y ))da(y), xeB d -\ n=l,2,- 
where 



fe=0 



Then V„/ £ nf„, KJ = / for all / £ n£, and for / £ ^(S^ 1 ) or C^- 1 ), 

(2.6) UK/ - /||p < cE n (f) p , l<p<oo. 

2.2. A class of differential operators on § d_1 . One of the main tools in our 
study is a class of differential operators Di.j, 1 < i j < d, that commute with the 
Laplace-Beltrami operator. Let SO(d) denote the group of rotations on R d , and 
let ei, ■ ■ ■ ,ed denote the standard orthogonal basis in R d . For 1 < i ^ j < d and 
i £ R, we denote by Qi,j,t the rotation by the angle t in the (xi, a^-plane, oriented 
such that the rotation from the vector ei to the vector ej is assumed to be positive. 
For example, the action of the rotation Qi,2,t £ SO(d) is given by 

(2.7) Qi,2,t(xi, ■ ■ ■ ,x d ) =(xi cost - x 2 sini, xi s'mt + x 2 cost, x 3 , ■ ■ ■ ,Xd) 

=(scos((?!) + t), ssin(0 + t),x 3 , ■ ■ ■ ,x d ), 

where (2:1,352) = s(cos</>, sin</>), and other Qi.^t are defined likewise. 

To each Q £ SO(d) corresponds an operator L(Q) in the space L 2 (S d_1 ), defined 
by L(Q)f(x) := f(Q~ 1 x) for x £ § d_1 , which is a group representation of SO(d). 
The infinitesimal operator of L(Qijj) has the form 

(2.8) A, 3 - I [L(g W|t) ] ^ = xj ± - x~, l<i<j<d. 

In particular, it is easy to verify that, taking — (1, 2) as an example, 

(0 \ r 
~d(j)) /( SCOS( ^' ssin< ? :) ' 2;3 '---' a;d )- 

The following useful observation, which asserts that D\ ■/ is independent of smooth 
extensions of /, is a simple consequence of (|2 .9[) : 

Proposition 2.1. Let £ and let F and G be two smooth functions on an 

open neighborhood U C M. d of xq which coincide on L/n§ d_1 . Then D^jF(xo) — 
DTjGixo). 

The operators Dij are connected to the usual tangential partial derivatives ac- 
cording to the following formula [3 ! , (3.15)]: 



< 2 ' 10 > 5J 



n r / \ n d 




= djf - Xj ^ X i d if = ~ X! X i D i,jf- 
IHI =1 i=l {i:l<4^j<d} 
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The operators Dij are also closely related to the Laplace-Beltrami operator Ao. In 
fact, Ao satisfies the following decomposition [3, (2.6)], 

(2.11) A = £ Dl y 

\<i<j<d 

Furthermore, each operator Dij in this decomposition commutes with Ao. In 
particular, by (|2.3I) . this implies that the spaces of spherical harmonics on S d_1 are 
invariant under D^j. 

2.3. Moduli of smoothness and K-functionals on S d_1 . For 1 < i < j < d 

and 9 G [— ir, ir], we define the rth difference operator A[^- g by 



k=0 



where Tq/(x) := f(Qx) for Q G SO(d). This differential operator can be expressed 
in terms of the usual forward difference as, for example, for = 1(, 2), 

(2.12) A r 12 e f(x) = Alf (x 1 cos(-) - x 2 sin(-), x 1 sin(-) + x 2 cos(-), x 3 , . . . , x d ) , 

where Ag is acted on the variable (•), and is evaluated at t = 0. The following new 
modulus of smoothness was recently introduced in [31 Definition 2.2]: 

Definition 2.2. For r e N, (> 0, and / 6 LP(S d ~ 1 ), 1 < p < oo, or f G C(S d-1 ) 
for p = oo, define 

(2.13) u r (f,t) P ■= sup max ||A[ - e /|| . 

|0|<ti<»<j<«*" p 

This modulus of smoothness enjoys most of the properties of classical moduli of 
smoothness ([3l Proposition 2.7]), and it permits both direct and inverse theorems 
(0 Theorem 3.4]): 

Theorem 2.3. For f G D>(S d - 1 ) if 1 < p < oo and f G C^ 1 ) ifp = oo, 

(2.14) E n (f) p < curfan,- 1 ^, l<p<oo. 
On the other hand, 

n 

(2.15) u: r {f,n- l ) p <cn- r ^k r - l E k _ 1 (f) p , 1 < p < oo. 

fc=i 

A new K-functional on S d_1 was defined in terms of the differential operators 
Dij in [3, Definition 2.4]: 

Definition 2.4. For r G N and i > 0, 



(2.16) K r (f,t) p := ^jrf 1 11/ - 9\\p + t r t max ^ ||D[ j3 || p 

As in the classical setting, these two gadgets are equivalent ([3J Theorem 3.6]). 

Theorem 2.5. Let r e N and Zet / G L^" 1 ) zf 1 < p < oo and f G C^ 1 ) if 
y> = oo. For < f < 1, 

Ur(f,t)p~K r (f,t) p , l<p<oo. 
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Finally, we point out that there are several well studied moduli of smoothness 
on (see, for instance, [8]). One of the advantages of our new modulus 

is that it reduces many problems in approximation on S d_1 to the corresponding 
problems of trigonometric approximation of one variable, the latter is classical and 
well studied. Another advantage is that our modulus is relatively easier to compute, 
as demonstrated in Part 3 of [3]. 

3. SOBOLEV SPACES AND SIMULTANEOUS APPROXIMATION ON S d_1 

The classical Sobolev space W£ on is defined via the fractional order 

Laplace-Beltrami operator (see, for example, [HOIS]): 

(3-1) W; := {/ G L^- 1 ) : \\f\\ w; := ||/|| p + || (-A r /2 /ll P < ^} ■ 

We shall introduce a new Sobolev type space on § d_1 in this section and then study 
approximation by polynomials for functions in this new space. Our new Sobolev 
space is defined via the differential operators -Dij, 1 < i < j < d, which, by (|2.1f p . 
are more primitive than Ao. First, however, we need a lemma. 

Lemma 3.1. For f,g G C^S** -1 ) and 1 < i / j < d, 

(3.2) f f(x)Di, jg (x)da(x) = - / D l , J f{x)g(x)da(x). 

Proof. By the rotation invariance of the Lebesgue measure da, we obtain, for any 

9 G [-7T,7T], 

/ f{x)g(Qi,j-ex)da(x)= f(Q iJy ox)g(x)da(x). 

Differentiating both sides of this identity with respect to 9 and evaluating the 
resulted equation at 9 — lead to the desired equation (13.21) . □ 

The equation (I3.2j) allows us to define distributional derivatives on § d_1 for 
r € N via the identity, 



D r l J(x)g(x)da(x) = (-l) r / f(x)Dl j9 (x)da(x), g £ C 00 ^ 1 ). 
s d -! Js d - 1 

We can now define our new Sobolev space on the sphere. 

Definition 3.2. For r £ N and 1 < p < oo, we define the Sobolev space = 
Wp(§ d_1 ) to be the space of functions f G L p (§> d ~ 1 ) whose distributional derivatives 
Dl jf , 1 < i < j < d, all belong to L p (S d_1 ), with norm 

\\f\\ w; ($*-i)--=\\f\\ P + E wv'iJWp, 

l<i<j<d 

where L p (E> d ~ 1 ) is replaced by C(S d ~ 1 ) when p = oo. 

The following proposition compares the new Sobolev space with the classical one 
defined in (l3~T1): 



Proposition 3.3. For 1 < p < oo and r — 1 or 2, one has 

W; = W r p and \\f\\ W r ~ \\f\\ w; . 
In general, for r > 3 and 1 < p < oo, 

W;cW r p and \\f\\ w; < \\f\\ W r. 
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Proof. These are immediate consequences of (3.13) and (3.17) of [3]. □ 
Theorem 3.4. IfreN, f e W p , and 1 < p < oo, then 
(3.3) E 2n (f) P <cn- r max E n {D\ ,/)„. 

l<i<j<d ,J 



Furthermore, V n f, defined by (|2.5p . provides the near best simultaneous approxi- 
mation for all -f ' , 1 < i < j < d, in the sense that 

(3.4) \\Dl } {f - V n f)\\ p < cE n (D r itj f) p , 1 < i < j < d. 

Proof. Applying (|3. 21) to the function g : [— 1, 1] h-> K and using 

D ij9{{^y)) = 9 , ((x,y}){x l y J - xjVi) = -D^g{{x,y)), 
it follows immediately that 

D\ J / f(y)9({x,y))da(y)= [ D id f(y)g{{x,y))da(y). 
Js d - 1 js^-i 

Consequently, by (f2~4|) and (|2~5j) , we see that V n D^ = DljV n . Thus, using Theo- 
rems 12.31 and 12.51 we obtain 

E-2nU)p = E 2 n(f - V n f) p < cK r (f - V n f, U'^p 

<cn- r max . ||£>r (f~V n f)\\ p 

l<j<]<d 

= cn~ r ma, -\\D\J- V n (D r t J) \\ p 

l<j<]<d 

< cn~ r m&x E n (Dl J) p , 

l<i<j<d J 

where we used (|2.16|) in the third step, the fact that V n D\ j = D\ jV„ in the fourth 
step, and (|2.6|) in the last step. This proves (I3.3[) . The inequality (|3.4[) follows 
immediately from the above proof. □ 

Next we define a Lipschitz space on the sphere and consider approximation in 
such a space. 

Definition 3.5. For r G N, 1 < p < oo, and a S [0,1), we define the Lipschitz 
space W r p a = W^ a {S, d - 1 ) to be the space of all functions f e W p with 

mm n-fii \W,j,e(Dijf)\\p 

where £ is a fixed positive integer, for example, I = 1 . 

Our next theorem gives an equivalent characterization of the space W p ' a . For 
the same set of parameters as in the definition of W p ,a , we define the space 

H r p +a ■= (/ e L^ 1 ) ■ \\f\\ H ; + * ■■= II/IIp + sup Ur+ ^ t)p < oo) . 
I p o<t<i i ^ J 

Theorem 3.6. IfreN, 1 < p < oo, and a e [0, 1), iften 



IKS 



POLYNOMIAL APPROXIMATION IN SOBOLEV SPACES 



7 



Proof. To prove that / G implies / G H r p +a and < c\\f\\ w ^, it 

suffices to show that for / 6 W p and £ EN, 

(3-5) ||A^/|| p <c|0n|A^. e (^./)|| p . 
Using Lemma 2.6 (ii) in [3], we have 

\\^!ef\\ P = \HoA A Lef)\\ P < c \ \ r \\ D U A Lef)\\ P - 

However, from (|2.12|) and ()2.9|) . a quick computation shows that Aij^Dij — 
DijAi^e, hence, by iteration, A^D^ = D^jA e ljg . As a result, 

Together, these two displayed equations yield (|3.5[) . 

Conversely, assume / G H r p +a . We first show that D r id f G L p (§ d " 1 ). For 
.g G C* 00 ^ 1 ) and G^(t) := g(Q i>jlt x), we can write [3 (4.8)] 



Jo J 
which implies, in particular, that 

(3.6) lim " J 't = GM(0) = Dl j9 {x). 



Thus, by the definition of the distributional derivative D\ ■/, it follow that, for 
g G C* 00 ^ 1 ), 



[^/(x)]ff(x) da = (-If / f(x)D r ij9 ( x ) da 



= (-iriim/ /(,)%^^ = (-iriim/ %^(,)^ 

where the last step uses the rotation invariance of the Lebesgue measure da. How- 
ever, by the Marchaud inequality (Proposition 2.7 of [3]), for / G Hp +a , 

Ur(f,t) p < C t f f Ur+ j f + f P du < cr\\f\\ H;+a . 



Hence, by Holder's inequality, we deduce with i + y = 1 that 



[D r l J(x)]g(x)da(x) 



<4f\U+49\\p> 



d-n 



which implies, upon taking supreme over all g with ||g|| P ' < 1 that D\ ■/ G L P (E>' 
Next we note that for I G N, 

(3-7) \KM d U)Wp < c f I! a ^/Hp4?t> 

Jo u 

which follows from the analogue result for trigonometric functions [H (7.1)] as in 
the proof of Lemma 2.6 of [3]. Consequently, it follows that 

||/|| wr <||/||p + c max sup f \\&*£j\\ p -£L < c\\f\\ H r +a 

l<i<J<d <|e|<l W\ Jo u 

since < a < 1. This completes the proof. □ 
Theorem 13.61 together with the Jackson theorem implies the following: 
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Corollary 3.7. If r G N, a G [0, 1), / G W£ a , andl<p<oo, then 

E n (f) p <cn- r - a \\f\\ w; ,«. 

4. Polynomial approximation on M d : Recent Progress 

In this section, we recall recent progress on polynomial approximation on M d as 
developed in [3J. For /j, > 0, let Wp denote the weight function on M d defined by 

(4.i) w^^a-wxfr- 1 / 2 . 

For 1 < p < oo we denote by ||/|| p ,p the norm for the weighted IP space L p (M d , W^), 

i/p 

(4-2) \\f\W-= [ I \f(x)\PW^(x)dx' 



and ||/||oo,/i := ||/||oo for / € C(M d ). When we need to emphasis that the norm is 
taken over B d , we write ||/|| Pi „ = ||/|U»> ( b^)- For / G LP{M d ,W^), 1 < p < oo, 
or / G C(E d ), p — oo, the best approximation by polynomials is defined by 

En(f) P ,n ■= inf. Wf-P\\w 

4.1. Weighted orthogonal polynomial expansions on M d . Let V d {W^) denote 
the space of orthogonal polynomials of degree n with respect to the weight function 
W p on B d . We denote by P£(x, y) the reproducing kernel of V*(W M ) in L 2 (M d , Wp). 
It is shown in 9, Theorem 2.6] that 

(4.3) P£(x,y) = j ^ i Z n . d+rn ((x,y) + Vl-\\y\\ 2 (x',0) dcrfc) 

for any x, y G M d and (x, x') G § d+m_1 , where Z n ^(€) is the zonal harmonic defined 
in (|2.4p and /i = ln ^~- For r\ being a C°°-function on [0, oo) that satisfies the 
properties as defined in Section 1.1, we define an operator 

(4.4) V r ?f(x):=a, [ f(y)KZ(x,y)W,(y)dy, x G B d , 



where is the normalization constant of Wp and 



2n f k\ 
K%(x,y) - 

fc=0 v J 



(4-5) K»{x,y):=Yn[-)PNx,y) 



This operator plays the same role as V n f in the study on S d x . In particular, 
V£f G U d n and ||/ - V£f\\ Ptll < cE n (f) p>fi , 1 < p < oo. 

The spaces V^(Wp) of orthogonal polynomials are also the eigenspaces of the 
following second order differential operator: 

d d 

(4.6) Dp := ^(1 - xt)d 2 - 2 £ x^c^ - (d + 2 M ) ^ x^. 

i=l \<i<j<d i=l 

Indeed, elements of V^(Wp) satisfy (cf. 6, p. 38]) 

(4.7) V tl P = -n{n + d + 2fi- l)P for all P G V^(Wp). 
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It was shown in [3j Proposition 7.1] that the differential operator can be de- 
composed as a sum of second order differential operators: 

(4-8) 2>, = I>M+ E D h= E D lv 

i=l l<i<j<d l<i<3<d 

where the operators Df ■ for 1 < i < j < d are defined as in (|2.8j) . and 

(4.9) £>?, ee Dl^ := [^(x)]- 1 ^ [(1 - Ikf)^^)] ft, 1 < * < d. 

For the rest of the paper, we will always set <p(x) = d\ — ||a;|| 2 . We have the 
following useful estimates: 

Proposition 4.1. J/1 < p < ao, and g £ C 2 (M d ) then 

(4.10) IMp,p~ E H^Hp./*- 

l<*<j<d 

Furthermore, i/rsN, 1 < p < oo, and 17 G C 2r (B li ) i/ien 

(4.11) Cl 11^0^11^ < HU^II^ < c 2 [|vJ 2r af r s|| P , P + c 2 ||5|| P , M , 1 < * < d. 

Proof. (|4.10[) was proved in [3l Theorem 7.3]. In the case when r = 1, (14. lip was 
shown in [3j Theorem 7.4], and the proof there works equally well for r > 1. □ 

4.2. Moduli of smoothness and K-functionals. Two moduli of smoothness 
and their equivalent K-functionals on M d were introduced in [3] . Our first modulus 
on M d was defined via an extension / of a function / : M d — > R to B d+1 , defined by 

(4.12) f(x,x d+1 ):= f(x), (i,%i)£l W , x£B d . 
More precisely, it is defined as follows (0 Definition 5.3]). 

Definition 4.2. For r £ N, t> 0, and f £ L p (M d , W^), 1 < p < 00, or f £ C(B d ) 
for p — go, define 

(4.13) uj r (f,t) Pifl := sup <^ max ||A[ . e /|| if , (B d w ); 

|0|<t L 1 ^ l< ^^ d 

™£X ||Aj )d+lje 7|| i „( ) Bd + i jW(i _ 1/2 )| , 

w/iere /or m = 1, ||A£ d+lje /|| iJ >( B .i+i > w A< _i /a ) is replaced by \\A r l d+i e f\\ LP{sd) . 

In the case when /x = and to G N, we have established in j3j Theorem 5.5] 
the direct theorem, that is, the Jackson inequality 

(4.14) En(f) P ,» < cw r (/,n -1 ) P)/1 , l<p<oo, 
and the corresponding inverse theorem, 

n 

(4.15) ur(f, n- 1 )^ <cn- r J2 k^MfW, 

k=l 

in terms of this new modulus of smoothness. Moreover, it was also shown in [3j The- 
orem 5.8] that the modulus of smoothness oj r (f, t) p ,^ is equivalent to the following 
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K-functional: 

(4.16) K r (f,t)p tfl := mi \\\f - g\\ LP( B^ w ^ 

+ f 1< 1 ^ d \\Di tj g\\ L p(B^w^ + t r max \\D r l d+1 g\\ LP[Kd+Kw ^ i/2) y 



where if m = 1, then \\D r i>d+l g\\ LV(fi d+i^ w ^_ y2) is replaced by ||-D[ >d+1 5|| L p( S *)- 

Our second new modulus in [3] can be considered as a higher-dimensional ana- 
logue of the classical Ditzian-Totik modulus on the interval [—1,1]. In the un- 
weighted case, this modulus is defined as follows [3j Definition 6.7]: 

Definition 4.3. For r £ N, t> 0, and f G L p (B d , W^), 1 < p < oo, or f £ C(B d ) 
for p — oo, define 

(4.17) u) r (f,t) p := sup \ max \\Al jh f\\ p max\\A r h f\\ 

0<\h\<t {l<t<J<d l<i<d 

where denotes the central difference and 

&h 9ei m ■■= E(-l) fe ([) /(* + (§- k)hp{x)ei) , 

and we assume A" l hipe . = if either of the points x ±r hv ^ does not belong to M d . 

The weighted version Q r (f,t) p ^ of the above modulus can also be defined, but 
is more complicated. Both direct and inverse theorems were established in terms of 
Q r (f,t) p in [31 Theorem 6.13]. Furthermore, it was shown in [31 Theorem 6.10] that 
the modulus of smoothness Q r (f,t) p is equivalent to the following K-functional 

K r (f,t) Ptli := inf - g\\ Ptfl + f max J-D? > . 5 || + f max : \\v> r dlg\\ p A 

geC r (M d ) v. l<i<j<d ' J l<i<d ) 

in the sense that, for the equivalence between K r {f,t) p and w r (/, t) p in the un- 
weighted case, 

c- 1 2 r (/,t) p <X r (/,t) p < C tD r (/,t) p + ct r ||/|| p . 

The two if-functionals, hence their equivalent moduli of smoothness, are con- 
nected as shown in [3j Theorem 6.2]. 

Theorem 4.4. Let fi = Zf± and m £ N. Let f £ L p (M d , W^) if 1 < p < oo, and 

f £ C(M d ) if p = oo. We further assume that r is odd when p = oo. Then 

(4-18) ^i(/,t)p,p~^i(/,t) P , M , 
and for r > 1, t/iere is a t r > suc/i £/iai 

(4.19) K r {f,t) Pifl <cK r {f,t) Ptll + cf\\f\\ Ptll , 0<t<t r . 

Finally, we point out that it was shown in [3J that both moduli D r (/, t) Pjft and 
w r(/, i)p,M en J°y m °st of the properties of classical moduli of smoothness and they 
are computable as demonstrated in Part 3 of _3. . In comparison, the only other 
modulus of smoothness [TO] on the unit ball that is strong enough to characterize 
the best approximation is hardly computable. 
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4.3. Representation of the term D\ d+1 f- The term D r id+1 g appears in the def- 
inition of our first if-functional K r (f, t) p ^ in (|4.16p on the ball, where g(x, Xd+i) = 
g(x) as in (|4.12p . Notice that / is a function in a; € R d , but the operator Dt^+i = 
Xidd+i — Xd+\di involves Xd+i, so that D\ d+1 f is indeed a function of (x, Xd+i) in 
jjd+i rp^ e f ]] owm g lemma gives an explicit formula of this term in terms of /. 

Lemma 4.5. Assume that (y,y d +i) = s(x,x d +i) G M d+1 with s = \\(y,yd+i)\\ > 0, 
x e M d and x d+ i = (p(x) > 0. If f E C r {M d ), then 



(Dld+J)(y, y d+ x) = { ) fax) 



,. , , 1 < i < d. 

Proof. The proof uses induction. For r — 1, we have 

Di,d+if(y,yd+i) = {Vidd+i -Vd+idi)f{y) = -y d +idif(y). 
Hence, using the fact that -^-[f(sx)} = s(dif)(sx) we have 

d 

(Di,d+if)(sx,sx d +i) = -sx d+ i(dif)(sx) = -sip(x)(dif)(sx) = -<p(x) — [f(sx) . 

OXi 

Let F r (x, Xd+i) = d+1 f(x, Xd+i)- Assume that the result has been established 
for r. Then F r (sx, Sif(x)) — (—(pdi) r [f(sx)]. By definition, 

(4.20) F r+ i(sx, sxd+i) = (D itd+ iF r )(sx,sx d+ i) 

= sXi(d d +iF r )(sx, sxd+i) - sx d+1 (diF r )(sx, sx d+ i). 

On the other hand, taking derivative by chain rule shows that 

d \ r+1 / d \ 

~' P ^dx : J ^( sx ^ = \ ip ^dx~- J i F r(sx, s(p(x))] 

= -sip(x)(diF r )(sx,sip(x)) + sxi(d d +iF r )(sx, sip(x)), 
which is the same as the right hand side of (|4.20j) with Xd+i = f{x). □ 

Lemma 4.6. The function F)\ d + if( x > x d+i) * s even in Xd+i if r is even, and odd 
in Xd+\ if T is odd. 

Proof. For r = 1, D iyd+ if(x,x d+ i) = -x d+1 dif(x) is clearly odd in x d+1 . And 

Di >d+ J{x,x d+ i) = -Xidif(x) + x d+1 dff(x) 
is even in x d +i- The general case follows from induction upon using (I4.20p . □ 

Recall that our if- functional K r (f,t) Pyll with /i = in (|4.16l) is defined, 

when m = 1, with \\Dl d+1 g\\ LP(§d) in P^ce of \\D r i d+1 g\\ LP(Bd +i tW ^_ i/2) . Hence, as 
a consequence of the above lemmas, we conclude the following: 

Proposition 4.7. For g 6 C r (M d ) and the Chebyshev weight Wo on B , we have 
( 4 - 21 ) \\F>ld+i9\\Lp^) = \\{^di) r g\\ LP{v ,^ Wo) . 
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Proof. Let §| = {i £ : Xd+i > 0}. By Lemma T4. 61 we only need to consider 
when dealing with d+1 g. By Lemma 14.51 with s = 1, we then obtain 

\ D l d+1 g(x,Xd+i)\ P dcr(x, Xd+i) = 2 / \(ip(x)di) r g(x)\ p da(x, x d+1 ) 

J&* 



which is what we want to prove. □ 
In general, using polar coordinates, Lemma l4.5[ and Lemma l4.6[ we can deduce 
Proposition 4.8. If g e C r (M d ), \i = and m > 1, then for 1 < p < oo, 
(4.22) \\Dl d+1 g\\ P LP{M ^^_ 1/2) 



s 



"(l-s 2 )^ 1 I \(<p(x)d i r\g(ax)]\ p u dX n da; 



whereas for p — oo, we have 



ma , x , \ D i,d+i9(v)\ = max 

x£B d ,0<s<l 



d V 
VWfaT. ) iaisx)} 

5. Sobolev Spaces and Simultaneous Approximation on B d 

We start with the definition of a Sobolev space on B d . 

Definition 5.1. For 1 <p < oo, f G C r (M d ), and r gN, we define 

d 

(5-1) ll/llw;(B W := II/IIp,/. + E \\ D iA» + E 11^/11^. 

l<i<j<rf i=l 

and define W£(B , W^) fee £/ie completion ofC r (M d ) with respect to the norm 
II ' llw;(B d ,iy M )- 

Remark 5.1. Since convergence in the norm || • ||w(B d ,H^) implies convergence in 
the weighted L^-norm || • || P)/1 , we may assume that W£(B d ,W M ) C L p (B d , ) 
when p < oo, and Wp(B d , W^) C C(B d ) when p = oo. As a consequence, we can 
also extend the definitions of the operators D\ •/ and f r dl f to the whole space 
W£(B d ,W^). 

The following proposition follows readily from (6.15) and (6.16) of [3] and Propo- 
sition PO 



Proposition 5.2. /// e C r (M d ), \i > and 1 < p < oo, ffeen 
( 5 -2) E ll- c *<,d+i/lli*(B d + 1 ,w^_ i ) < c \\f\\w;(s d ,w^)- 



Ki<d 



Furthermore, if f & C (B ) and 1 < p < oo then 
(5-3) ll/llw=»-(B<«,Wi,) ~ E \\ D Zjf\\p^- 

l<i<i<d 
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Theorem 5.3. Let fi = **=± mth m G N. For f G W;(B d , W^), 1 < p < oc, 



(5.4) E 2n (f) P!fl < 



cn 



max E n {D\ jf)^ 
' P} B ^ J E n(Dl d+1 f) LP md + i w 

Ki<d ' 



Furthermore, V^f, defined by (|4.4[) , provides the near best simultaneous approxi- 
mation for all D\ jf , 1 < i < j < d + 1 in the sense that 

\\DUf - V£/)|| Pj „ < cEniDTJ)^ l<i<j<d 

\\Di,d+l(f ~ ^"/)||.LP(B d + 1 ,W _i) < c£; n(A I '.<i+l/)LP(B t! + 1 ,Vl/ i)j 1 < i < d. 

Proof. For / defined on B d , we define F(x,x') := f(x), x G M d , (x,x') G g^™- 1 . 
By |1 Lemma 5.2], (V^F)(i, x') = V£f(x). Furthermore, by Lemma 5.7], 
K r (f,n^ 1 ) LP ( K d W ^ ~ K r (F, n~ 1 ) LP ( S d+ m -iy Hence, it follows that 

£ , 2n(/)iP(B< 1 ) < Ci^ r (/ - l^f/, " _1 )LP(B d ,W M ) 

< cK r (F - V„f, n _1 ) iP ( St i+ m -i) 

< cn~ r max 1 1 D\ , (P - F„F) || iJ>(sd+ m-i 3 



cn^ max ||£>[ -F - (£»[ ,-V„F|| L p (s h 

Ki<j<d+1 



where the last step follows from the fact that V n (D^ d+k F) depends on Xj, 1 < j < d, 
and Xd+k , which implies that we only need to consider l<i<j<d + l. 

Denote by Vff d the operator (|4.4|) associated with on B d and f(x, xa+i) = 
f(x). By Lemma 5.2 of [3], 

(5.5) V n F(x, x') = V»- d l{ 2 J{x, x d+l ) = Vl J(x). 

Since D\ -V n — V n D"[ ■ on the sphere, it follows that, for 1 < i,j < d, 

DljV^fix) = D\j(y n F)(x,x') = V n D r i j F(x, x') = V^DVJix). 
Consequently, it follows from [3j (5.8)] that 

\\D r itj F - DjjV n F)\\ LP (§d+m.-i) = 4DVJ - DIXJK^Wlhm^w.) 

= 4Dijf-V£ i (P r id f)(x)\\ L ,p*, Wll) 

< cE n (Dljf) L p^d, Wll y 
Whereas for D\ d+1 F term, we have for 1 < i < d, 

Dld+xVZdlif^d+x) = Dl d+1 (V n F)(x,x') 

= V n (Dl d+1 F)(x,x') = V^DlJ^Xd+x). 

Consequently, 

\\ D l,d+i F ~ D l d+1 V n F\\ LP ( S d +m -i) = c\\D r itj f - V^ d l{ 2 Dl j f(x)\\ LP(B d+i. W)i _ i/2 ) 

This proves (|5.4j) . The conclusion that V£f is the near best simultaneous approxi- 
mation follows from the above proof and (|5 .5[) . □ 
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Corollary 5.4. Let p = Z*=± with m G N. /// G W£(B d , W M ), 1 < p < oo, i/ien 

E n (f)p,n < cn~ r \\f\\ W r(Md,w„)- 

Proof. This follows immediately from (|5.1I) . (|5.2p and Theorem 15.31 □ 

In the next corollary, we replace D\ d+1 f term in (|5.4j) by ordinary derivatives 
of /. First we consider the Chebyshev weig ht Wo on B d (with p, = Q). 

Corollary 5.5. For f G W£(B d , VFo), r G N and 1 < p < oo, 

E2n(f) P .o < cn~ r max E n (D\ J) p + cn _r max {{(fidi) r f) . 

l<l<j<a l<i<a r ' 

Froo/. By (jl^Tj) . for all / G C r (B d ), 

(<^(a;)a 4 ) r /(a;) = D[ id+x /(a;,a;d + i), 

where x GM d and a^+i = vK 3 -)- The desired conclusion then follows. □ 

For /i > 0, including the case p = 1/2 (the constant weight function), however, 
the best that we can do is the following: 

Corollary 5.6. Let p, = 2^1 and rn G N. For / G W;(B d , 1 < p < oo, 
E 2n {f)p,ti < crT r max E n {D\ ,•/)«,« + cn~ r max max £7„_ r (Q?/)„ jP 



max E n -r(^/)p, M +0/-§)p 

^^-<j<r 



Proof. It was shown in Lemma 6.4 of [3] that 



G B d+1 , 



3=1 

where pj !r is a polynomial of degree < j. Since E n (f) is subadditive, it follows that 

r 

(5.6) E n (Dr d+1 /) iP(B(i+liW(i _ i/2) < ^ inf PjA d if-g) 



L»(l i + 1 ,W„_ 1/2 ) 

However, using (6.11) and (6.12) of [3J, we have, for (a;, x^+i) G B d+1 , 



|pj, r (a;,a; d+ i)| < 



c, ifl<j<2±i, 
ckd+i| 2j - r , if^<j<r. 



Thus, by (|5.6p . we deduce 

- B n(-D[,d+i/)Lf(B d + 1 ,w M _ 1/2 ) < c max inf \\dff-g\\ p ^ 

+ c max inf / - g\\ p ^ +u _- )p . 

The desired conclusion then follows from Theorem 15.31 □ 



It remains to be seen if D\ d+1 f term in (|5.4p can be bounded by a term that 
involves only (ipd) r f in the case of p, > 0. 

Similar to the case of we can also define a Lipschitz space on the ball. 
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Definition 5.7. For r £ N, a £ [0, 1),. and 1 < p < oo, we define W£ a (M d , W M ) to 
be the space of all functions f : M d — > K wii/i finite norm 

\\f\\w;>"(m W ■= + i max sup |^1A^(Z>[ /)|| P)/J 

+ max sup |6»r Q ||Af d+1 , 9 (-D[ d+1 /)|| Lt , (B<J w _ 1/2 ) 

l< l <d 0< | 9 |< 1 

wif/i the usual change when p = oo, where £ is a fixed positive integer, say £ = 1. 

We can also give an equivalent characterization of the space Wp a (M d , W^) in 
terms of our modulus of smoothness. For the same set of parameters as in the 
definition of Wp ,a (B d , we define a space 

Hl+«(M d ,W,) := (/ G LP(B d , W») : sup < oo) . 

I o<t<i i J 

Theorem 5.8. Let fi = ^f^. IfreN, 1 < p < oo, and a £ [0, 1), then 

W;' Q (B d ,Wy = Hp +r (M d , W M ) and ||/|| w ;.« (Bd) ^) ~ ||/||^ (b ,, Wm) . 

Proof. This follows from Theorem EH since for F(x,x') := f(x), (x,x') £ S^" 1 " 1 
and x £ B d , we have u>(F, i)j>(§<i+m-i) ~ uj r {f 1 t) P4L by Lemma 5.4 of Ej. □ 

We could also define a Lipschitz space that uses central differences of &[ f in place 
of D ld+J in W£' Q (B d , WJ, so that it is equivalent to an analogue of H 7 p +a (M d , W^) 
with bJ r+ i(f, in place of uj r +i{f, t) p ^. 

As a consequence of the last theorem and the Jackson estimate, we have 

Corollary 5.9. Let fi = and m £ N. If r £ N, a £ [0, 1), / G W£ a (R d , Wp), 
and 1 < p < oo, then 

(5-7) E n {f) p ^ < cn- r - a \\f\\ w r, a(mitWfi) . 

Let us point out that for / G C r (M d ) the traditional definition of the Lipschitz 
continuity takes the form, for < a < 1, 

(5.8) \d p f(x)-d p f{y)\<c\\x-y\\ a , (3 £ N d , |/3| = r 

for all x, y £ M d . Let us denote by Lip r a the space of all C r (M d ) functions that 
satisfy (|5.8j) . From the definition of Dij it follows readily that 

Li Pr>Q c 

Hence, the estimate ()5.7j) holds for the functions in Lip rQ . On the other hand, our 
definition of Wp' a is more general than Lip r a as the following example shows. 

Example. Let f a (x) = (1 - ||x|| 2 + \\x - x \\ 2 ) a on M d with a fixed x £ S d_1 . 
Assume 1/2 < a < 1. Then by [3j Example 10.1], u) r (f a ,t)oo ~ * 2q , so that by 
Theorem 5.8, / Q G W^ 2a - 1 (B d ). ' On the other hand, setting x = (1, 0, 0, . . . , 0) 
shows that / Q (a;) = (1 — x\ + (1 — xi) 2 ) a = 2"(1 — xi) a , whose first partial derivative 
is unbounded on M d so that it is not an Lip 1 a function. We note that Di,jf a £ 
C{M d+1 ) for all 1 < i < j < d+ 1. 
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